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The internal structure of the thin transonic disk accreting onto a nonrotating 
black hole inside the last stable orbit (r < 3rg) is considered within the hydrody- 
namical version of the Grad-Shafranov equation. It is shown that in the vicinity 
of the sonic surface takes place a sharp diminishing of the disk thickness. As a 
result, in the vertical balance equation the dynamical forces p[(vV)v]g become 
important, these on the sonic surface being of the same order as the pressure 
gradient Ve/P. In the supersonic region the thickness of the disk is determined 
by the form of ballistic trajectories rather than by the pressure gradient. 


1 Introduction 

According to the standard disk model [jl| , the matter forms a thin balanced disk 
and performs a circular motion with keplerian velocity fi<;(r) = (GAt/r) 1 ^ 2 . 
The disk is thin provided that the accreting gas temperature is sufficiently low 
(c s -C i’k), so that H ~ rc s /v k- Introducing the viscosity parameter agg < 1, 
relating the stress tensor t r v and the pressure as t r v = «ss-P [0 , one can obtain 

v r /v K « a ss c 2 Jv^. (1) 

General relativity effects result in two important properties: the absence 
of stable circular orbits at small radii (r < 3r g for a nonrotating black hole, 
r g = 2 GM./c 2 is the gravitational radius), and the transonic regime of accre¬ 
tion. The former means that for r < 3r g a flow can be realized in the absence 
of viscosity. The latter results from the fact that according to ([[]) the flow 
is subsonic outside the marginally stable orbit r = 3r g while at the horizon 
r = r g the flow is to be supersonic. 

Up to now in the majority of works the procedure of vertical averaging 
was used, where vertical four-velocity u§ was assumed to be small |3J . As a 
result, the vertical component of the dynamic force nu b Wb(pu§) was postulated 
to be unimportant up to horizon. Here we are going to demonstrate that the 
assumption = 0 is not correct. As in the Bondi accretion, the dynamic 
force is to be important in the vicinity of the sonic surface. 


2 Basic equations 


Consider ideal gas accreting onto a black hole inside the marginally stable 
orbit r = 3r g . Exact equations of motion of ideal media in Kerr metric were 
formulated in (4|]. We consider a non-rotating (Schwarzschild) black hole only. 

In Boyer-Lindquist coordinates (t, r, 9, ip) the Schwarzschild metric is [|l| 
ds 2 = — a 2 dt 2 + dx l dx k , where a 2 = 1 — ‘IM.fr , g rr = a~ 2 , ggg = r 2 , and 
g w = w 2 = r 2 sin 2 9. Here indices without caps denote vector components 
with respect to the coordinate basis d/dr, djd9, d/dp and indices with caps 
denote their physical components. The V& symbol always represents a covari¬ 
ant derivative in space with g^ metric. Finally, below we always use a system 
of units with c = G = 1. 

ft is convenient to introduce a stream function 4>(r, 9). This function defines 
the physical poloidal four-velocity component u p as 

cmup = —(V4> x e^). (2) 
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Here n is the particle concentration in the comoving reference frame, ft is the 
curves $(r, 9) = const that define lines of flow of the matter. 

Three integrals of motion are conserved on $(r, 9) = const surfaces: en¬ 
tropy S = 5'($), energy F(<f>), and ^-component of angular momentum L(4>) 


F(<f>) = pay, L(<1>) = gwup. (3) 

Here /i = ( p m + P)/n ( p m is internal energy density) is the relativistic enthalpy. 
Below for simplicity we use the polytropic equation of state P = k(S)n r , so 
that the temperature and the sound velocity can be written as [!lj 

T = k(S)n T -\ c 2 = -k^S/n^ 1 . (4) 
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As a result, the relativistic Euler equation 


nu b V b (yiiua) + V a P - gn(up) 2 —V a w + nnj 2 —\7 a a = 0 (5) 

w a 

can be rewritten in the form of the Grad-Shafranov scalar equation for the 
stream function 4>(r, 9) containing three integrals £’($), T(<f)), and 5'(<f > ) |I| 
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in which the V' k derivative acts on all the variables except the quantity M 2 . 
Here the thermodynamical function M 2 is defined as M 2 = 47r/z/n and 
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To make the system complete we need to supply Grad-Shafranov equation (]6|) 
with the relativistic Bernoulli equation u 2 = q 2 — — 1; the latter with the 

help of (HI) can be rewritten as follows: 


2 (m2 2 r 2 / 2 2 2\ / 2 2 

u p = (E — a L jm — a fi )/ot /j . 


( 8 ) 


3 Subsonic flow 


First of all, let us consider the subsonic region in the very vicinity of the 
marginally stable orbit r = r 0 = 3r g where the poloidal velocity is much 
smaller than that of sound. Then equation (j6|) can be significantly simplified 
by neglecting the terms proportional to D ~ u 2 /c 2 . As a result, we have 
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This equation describing the subsonic flow is elliptical. Hence, it is necessary to 
specify five boundary conditions on the surface of the last stable orbit r 0 = 3r g 
where a 0 = a(r 0 ) = x^/3, u^(r 0 ) = l/v 7 ^, and y 0 = j(r 0 ) = \/4/3 [0. We 
consider below the case where the radial velocity is constant on the surface 
r = r 0 and the toroidal four-velocity is exactly equal to u^,(r 0 ): 

u f (r o, 0) = -M 0 , Mg(ro, 0) = 0u o , w^(r 0 , 0) = l/Vs. (10) 

Here Uq (r 0 , 0) corresponds to the plane flow at the marginally stable orbit, 
and we introduced the new angular variable 0 = 7 t/2 — 8 (©disk ~ c 0 ) which is 
counted off from the equator in the vertical direction. Next, we suppose that 
the velocity of sound is also a constant on the surface r = r 0 : 


c s (H),0)=Co. (11) 

For P = k(S)n r this means that both the temperature T 0 = T(tq) and the 
relativistic enthalpy /r 0 = /i(r 0 ) are also constant on this surface. It is necessary 
to stress that, according to (0), for nonrelativistic temperature of the accreting 








gas c 0 < 1 we have a small parameter u 0 /c 0 ~ ass c o •C 1. Finally, as the last, 
fifth, boundary condition it is convenient to specify the entropy S(Q). 

Introducing now the values eo = «o7o = \/8/9 and lo = M^(r 0 )r 0 = \/3r g , 
one can rewrite the invariants E(<&) and L(<E>) as 

E($) = /i 0 eo = const, L($) = fj, 0 l 0 cos Q m . (12) 


Here 0 m = 0 m ($) is the angle for which $(r 0 , 0 m ) = $(r, 0). In other words, 
the function Q m (r, 0) has the meaning of a theta angle on the last stable 
orbit connected with a given point (r,@) by a line of flow <f>(r, 0) = const. In 
particular, 0 m (r o , 0) = 0. 

First of all, we see that condition E = const ( 0 ) allows us to rewrite 
equation flf]) in a simpler form 


<9 2< 3> cos© <9/1 <9$\ 
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-W^^Wcos 2 ©^. (13) 
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Next, as one can show §, for r = r 0 , the r.li.s. of equation (p~3|) describes the 
transverse balance of a pressure gradient and a gravitational force, whereas 
the l.h.s. corresponds to the dynamic term (vV)v. At the marginally stable 
orbit it is of the order of Uq/cq and may be dropped. It is therefore natural 
to choose the entropy S(<$>) from the condition of a transverse balance on the 
surface r = r 0 
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where L(0 TO ) is determined from the boundary condition (|T^). Thus, we have 
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Owing to flU]), relation flL5|) corresponds to the standard concentration profile 


n(r 0 , 0) « TT-o exp 1 • ( 16 ) 

Finally, definition (]2j) results in the following relationship between functions 
$ and 0 m : 

d<E> = 27ra 0 ron(r 0 , <d m )u 0 cos 0 m d0 m . (17) 

Hence, due to (|T2|), ([Dd]), and (0), the invariant L(<h) can be directly deter¬ 
mined from the boundary conditions as well. 








Equation (0) together with boundary conditions (0), (ig), ([3, and 
([17]) determines structure of the inviscid subsonic flow inside the marginally 
stable orbit. For example, for qCl we obtain using (0) 
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depending on the radius r only, is a poloidal four-velocity of a free particle 
having zero poloidal velocity for r = rq. As we see, w 2 increases very slowly 
when moving away from the last stable orbit. Therefore, the contribution of 
w 2 turns out to be negligibly small in the subsonic region. 

An important conclusion can be drawn directly from (p~8D in which for the 
equatorial plane we have 0 m = 0 = 0. Assuming u p = c s = c* and neglecting 
w 2 , we find the velocity of sound c* on the sonic surface r = r*, 0 = 0: 


c* ~ \j C °' ( 20 ) 

As we see, c* ~ cp. Next, as the entropy S remains constant along the flow 
lines, the gas concentration remains approximately constant along the flow 
lines (n(r*,@) ~ n(ro,0 ra )) as well. In other words, in agreement with the 
Bondi accretion, the subsonic flow can be considered incompressible. On the 
other hand, because the radial velocity increases from up to c* ~ Cp, i.e., for 
Co < 1 ( Up/cp -C 1) it changes over several orders of magnitude, the disk 
thickness H should change in the same proportion owing to the continuity 
equation (see Fig. 1]) 

H(r *) « -H(r 0 ). (21) 

Co 

As a result, a rapid decrease of the disk thickness should be accompanied by 
the appearance of vertical component of velocity which also should be taken 
into account in Euler equation (|5|). 

Indeed, as one can End analyzing asymptotic of equation (|l3|) ||, in the 
vicinity of the sonic surface located at r* = rp — Au^ 3 rp, where the logarithmic 
factor A = (3/2) 2//3 [ln(co/no)] 2 ^ 3 ~ 5 — 7, the components of the velocity and 
the pressure gradient in the limit r —> r* can be presented as 










Figure 1: The structure of a thin accretion disk (actual scale) for c 0 = 10 -2 , 
Uq = 1CT 5 after passing the marginally stable orbit r = 3 r g . In the vicinity of 
the sonic surface r = r* the flow has a form of an ordinary nozzle. 


On the other hand, near the sonic surface the radial scale 5r determining 
the radial derivatives becomes as small as the transverse dimension of a disk: 
5r ps if(r*) ~ Uor 0 . Hence, logarithmic derivative rj\ = ( r/n){dn/dr ) can be 
evaluated as rfr ps Uq 1 . As a result, both components of the dynamic force 

Uq Ouq eg O . Cq Q r9 ox 

r 90 u 2 r ’ r 9r u 2 r ’ 1 J 

do become of the order of the pressure gradient. 


4 Transonic flow 


To check our conclusion one can consider flow structure in the vicinity of the 
sonic surface in more detail. Because the smooth transonic flow is analytical 
at a singular point, one can write 


n = n* ^1 + ijih+ ^h 3 © 2 + • • -J , (24) 

@ m = CLq f 0 + CL± h® + —Cl2h 2 0 + —6 q0 3 + • • • 


( 25 ) 


























where h — (r — r*)/r*. Here we assume that all the three invariants E, L, and 
S are already given. Hence, the problem needs one extra boundary condition. 
Now comparing the appropriate coefficients in Bernoulli (| 8 ]) and full stream 
equation (| 6 |), one can obtain neglecting terms ~ «o/ c o 
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where a 2 = ct 2 (r*) rs 2/3. 

As we see, coefficients (0)-(|3(]) are expressed through the radial loga¬ 
rithmic derivative rj\. They have clear physical meaning. So, gives the 
compression of flow lines: a o = H(ro)/H(r*). In agreement with ([H|) we 
have a 0 ~ c 0 /u 0 . Further, aq corresponds to the slope of the flow lines with 
respect to the equatorial plane. As aq > 0 , the compression of stream line 
finishes somewhere before the sonic surface, so inside the sonic radius r < r* 
the stream lines diverge. On the other hand, as cq -C Uq , for r = r* the di¬ 
vergency is still very weak. Hence, in the vicinity of the sonic surface the flow 
has a form of an ordinary nozzle (see Fig. 1). Finally, as 02 ~ i ]3 ~ bo ~ Uq 2 , 
one can conclude that the transverse scale of the transonic region //(r*) is 
the same as the longitudinal one. The latter point suggests a very important 
consequence that the transonic region is essentially two-dimensional, and it is 
impossible to analyze it within the standard one-dimensional approximation. 


Let us stress that it is rather difficult to connect the sonic characteristics 
hi = hi( r *) with physical boundary conditions on the marginally stable orbit 
r = r 0 (for this it is necessary to know all the expansion coefficients in ([?|) 
and (^5D). In particular, it is impossible to formulate the restriction on five 
boundary conditions (0), (0) and (|H>|) resulting from the critical condition 


on the sonic surface. Nevertheless, the estimate rft ~ Uq 1 makes us sure that 
we know the parameter r/i to a high enough accuracy. Then, according to 
(|30|), all the other coefficients can be determined exactly. 

Using now expansions (24) and (25), one can obtain all other physical 








parameters of the transonic flow. In particular, we have 
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Hence, a shape of the sonic surface u p = c s has the standard parabolic form 


h = 



(31) 


5 Supersonic flow 


Because the pressure gradient becomes insignificant in supersonic region, the 
matter moves here along the trajectories of free particles. Neglecting VgP 
term in the d-component of relativistic Euler equation (^|), we have (cf. 0) 
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Here can be easily expressed in terms of radius: = V3/x, where x = r/r g . 

We also introduce dimensionless functions f(x) and g(x): Qf(x ) = xUq, 
g(x) = —auf > 0. Using now (^) and definitions above, we obtain an ordinary 
differential equation for /(x) 


d / / 2 + 3 
dx x 2 g(x) 


(33) 


Next, analyzing equation (|D, one can conclude that u p — > w as r — * r g . On the 
other hand, u p ~ c* ~ Cq for r < r*. Therefore, the following approximation 
should be valid throughout r g < r < r* region: g(x) ~ \J(aw) 2 + (ac 0 ) 2 , 
where, owing to (|I9|), (ccw) 2 = (3 — x) 3 /(9x 3 ). 


The results of calculations are presented on Fig. |l[. As we see, in the 
supersonic region the flow performs transversal oscillation about the equatorial 
plane, their frequency independent of the amplitude. Once diverged, the flow 
converges once again at the point r = r cr (r 0 — r cr = Ac*ro, A ps 27t) on the 
equatorial plane. Such oscillation can be easily understood. Indeed, as one can 
see from (IB), for Co 1 poloidal motion remains nonrelativistic for r ~ r*. 
Hence, it is possible to use the nonrelativistic equation 
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where Vq = rdQ/dt. Because the spatial amplitude of oscillations ~ c 0 r 0 
is small compared with r 0 , one can write v$ = u ${ tq ) / 70 = 1/2 (see O), 
r = 3 r g , and |ty| ~ c*. This gives ©(h) = ©a sin[(n^/c*)h], h = (r — r*)/r* 
in full agreement with exact calculations. Clearly, additional consideration is 
necessary to determine the flow structure for r < r cr . Nevertheless, one can 
be sure that the accretion disk thickness, though oscillating in the supersonic 
region, remains as small as in the region of stable orbits: ©disk 1$ c 0 . 

6 Conclusion 

As was shown, the diminishing disk thickness in the vicinity of the sonic surface 
inevitably leads to the emergence of the vertical velocity component of the 
accreting matter. As a result, the dynamic term (vV)v in the vertical balance 
equation cannot be omitted. It is necessary to stress that whereas at the 
sonic surface both components of dynamic term [(vV)v]j (^) become of the 
same order of magnitude as the pressure gradient, the role of the gravitational 
term remains unimportant: Vg(p g ~ @/r, i.e., it is Cq/uq times smaller than the 
leading terms. As a result, the structure of a thin transonic disk is quite similar 
to that of an ordinary planar nozzle. For this reason the critical condition on 
the sonic surface does not restrict the accretion rate. We suppose that for a 
given accretion rate it determines vertical component of the velocity on the 
marginally stable orbit which does not affect the flow structure. Finally, inside 
the sonic surface r < r* the pressure term V gP becomes unimportant, so the 
thickness of the disk is determined by the form of ballistic trajectories. 

This work was partially supported by the Russian Foundation for Basic 
Research (grant No. 00-15-96594). 


[1] S.L. Shapiro and S.A. Teukolsky, Black Holes, White Dwarfs, and Neutron 
Stars, A Wiley-Interscience Publication, New York, 1983. 

[2] B. Paczynski and G.S. Bisnovatyi-Kogan, Acta Astron. 31 (1981) 283. 

[3] M.A. Abramowicz, A. Lanza, and M.J. Percival, Astrophys. J. 479 (1997) 
179. 

[4] V.S. Beskin and V.I. Pariev, Physics Uspekhi 36 (1993) 529. 

[5] V.S. Beskin, R.Yu. Kompaneetz, and A.D. Tchekhovskoy, Astron. Lett. 28 
(2002) 543. 


